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We report inelastic neutron scattering measurements of crystal field transitions in PrFeAsO, 
PrFeAsO0.87F0.13, and NdFeAsO0.85F0.15. Doping with fluorine produces additional crystal field 
excitations, providing evidence that there are two distinct charge environments around the rare 
earth ions, with probabilities that are consistent with a random distribution of dopants on the 
oxygen sites. The 4/ electrons of the Pr 3+ and Nd 3+ ions have non-magnetic and magnetic ground 
states, respectively, indicating that the enhancement of T c compared to LaFeAsOi-^Fj, is not due 
to rare earth magnetism. 



Tuning materials to enhance their properties is the 
driving force behind much of modern condensed mat- 
ter physics. Chemical doping is often the most prac- 
tical means of accomplishing this goal but a thorough 
microscopic understanding of the effect of doping is a 
major challenge. For example, the role of phase separa- 
tion in the high-temperature superconducting cupratesp] 
and colossal magnetoresistive manganites 2J is still being 
debated. In the recently discovered iron-based supercon- 
ductors, chemical doping is the primary method of induc- 
ing superconductivity, but there is conflicting evidence 
whether this is due to chemical pressure [3], the change 
in carrier concentration 4] , disorder [5], or a combination 
of all three. Moreover, the question of phase separation 
in the superconducting phase itself has not been con- 
clusively answered and may be material-dependent, with 
evidence of both phase separation into antiferromagnetic 
and superconducting regions [6] and phase coexistence [7]. 

The ^FeAsOi^Fz (R =La, Ce, Pr, Nd, Sm, Gd) 
series [SHIS] was the first family of iron-based supercon- 
ductors to be discovered. Replacement of lanthanum 
with other rare earths increases T c up to ~ 55 K in the 
optimally doped regime close to x = 0.15. Doping with 
fluorine adds electrons into iron ci-bands but the substi- 
tution of trivalent rare earths for lanthanum does not 
change the carrier concentration, so any changes are ei- 
ther due to the influence of local 4/ magnetism or the ef- 
fect of chemical pressure from the well-known lanthanide 
contraction. Understanding which of these are respon- 
sible for the enhancement of T c has not been resolved. 
Recent neutron diffraction, muon spin relaxation (/iSR), 
and Mossbauer studies of the i?FeAsO parent compounds 
have shown evidence for strong coupling between rare 
earth and iron magnetism [TMT^] . Assuming that this 
coupling persists in fluorine-doped systems, we might ex- 
pect the 4/ moments to influence the superconducting 
properties as well. Indeed, a recent 19 F NMR study of 



SmFeAsOi_ x F x has inferred a non-negligible coupling be- 
tween the 4/ and conduction electrons [20). 

Measurements of rare earth crystal field (CF) excita- 
tions using inelastic neutron scattering provide unique 
insight into the effects both of substituting magnetic rare 
earth ions for lanthanum and of doping fluorine onto the 
oxygen sublattice. This is because it is both a local probe, 
since the CF potential acting on the 4/ electrons is de- 
termined by the electrostatic environment produced pre- 
dominantly by the nearest-neighbor oxygen/fluorine ions, 
and a bulk probe, since the CF transition intensities are 
a true thermodynamic average of the whole sample. Rel- 
ative changes in the CF peak intensities can be directly 
related to the volume fraction of rare earth ions affected 
by a particular configuration of neighboring ligands. 

In this article, we present a comparison of CF excita- 
tions measured in PrFeAsOi-^F^with x = 0.0 and 0.13, 
which shows that there are two distinct charge environ- 
ments in the superconducting compound, similar to the 
conclusions of a recent 75 As NQR study However, 
the measured reduction in CF intensity with x is consis- 
tent with a random distribution of fluorine ions, show- 
ing that the charge environments are produced by the 
dopant ions and not electronic phase separation in the 
iron layers, as proposed in the earlier work. Similarly, our 
measurements on NdFeAsO0.85F0.15 reveal more crystal 
field excitations than would be allowed due to either 
tetragonal or orthorhombic symmetry, also indicating 
the presence of inequivalent rare earth sites. Moreover, 
we are able to eliminate or severely constrain the rele- 
vance of rare earth magnetism to the pairing mechanism 
since superconducting PrFeAsO0.s7F0.13 has a singlet 
ground state while superconducting NdFeAsO0.85F0.15 
has a magnetic Kramer's doublet ground state, even 
though both the superconducting transition tempera- 
tures are similar [51 110). 

Powder samples of i^FeAsOi^F^were synthesized fol- 
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lowing the method described in Ref. [T5]. Supercon- 
ducting transition temperatures determined by the onset 
of diamagnetism in an applied field of 20 Oe are 41 K 
and 49 K for PrFeAsOi^F* (x = 0.13 ±0.01 determined 
from the phase diagram of Ref. [H]), and NdFeAsOF .i5 
(nominal composition), respectively. Structural charac- 
terization by neutron diffraction was performed on GEM 
at the ISIS Facility and HB2A at the High Flux Isotope 
Reactor. Inelastic neutron scattering (INS) studies were 
conducted on time-of-fiight spectrometers Merlin at ISIS 
and IN4 and IN6 at the Institut Laue Langevin. The INS 
data have been placed on an absolute scale by normal- 
ization to a vanadium standard. 

Fig. 1 shows the neutron diffraction pattern for 
PrFeAsO0.87F0.13 measured at 4K, compared to a Ri- 
etveld fit using FullProf [23] • Structural refinements con- 
firmed the orthorhombic (Cmma) structure in the par- 
ent compound and the tetragonal (P4/nmm) structure in 
the superconducting sample, with evidence for 4.3±0.2% 
of FeAs impurity in both the x = 0.0 and 0.13 com- 
pounds. A PrOF impurity phase (~ 3.8%) is also de- 
tectable in PrFeAsO0.87F0.13, but is too small to pro- 
duce the CF excitations discussed below. There is no ev- 
idence of structural inhomogeneity in the primary phases 
of PrFeAsO and PrFeAsO0.s7F0.13. No impurity phase 
was observed in NdFeAsO0.85F0.15. 

Fig. 2 shows the CF excitations in PrFeAsO and 
PrFeAsO0.87F0.13 below 14meV. The crystal field exci- 
tations are evident as the peaks at nonzero energy trans- 
fers that are not present in LaFeAsO, which has no 4/ 
electrons, There are additional CF excitations between 
30 and 40 meV, which are not shown. In the par- 
ent compound, PrFeAsO, the CF excitations centered 
at ~ 3.5 meV are split above TN(Pr)=12K. We assume 
that this is due to the internal molecular field created 
by the Fe sublattice below T N (Fc) = 127KQ2], although 
we cannot confirm that they become degenerate above 
T]sr(Fe) because of thermal broadening. We have insuffi- 
cent information to solve the CF potential so we cannot 
construct a microscopic model to explain the collapse of 
the splitting when the rare earth sublattice magnetically 
orders as seen in the inset of Fig. 2(a). However, it may 
be due to a spin reorientation, similar to what has been 
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FIG. 1: Neutron powder diffraction in PrFeAsO0.87F0.13 mea- 
sured on HB2A at 4K (circles) compared to the Rietveld re- 
finement (line). The line below the plot shows the difference 
between the data and the refinement. 



observed in NdFeAsO[23], and is consistent with other 
evidence of an interplay between iron and rare earth mag- 
netism in PrFeAsO, such as a reduction in the intensity 
of the iron magnetic Bragg peak|16] and a reduction in 
the /xSR frequency [IB] approaching Tn(Pi"). 

In PrFeAsO0.87F0.13, there is also a CF peak at ~ 
3 meV (labeled A) but there is no evidence of any split- 
ting, which is consistent with the absence of long range 
magnetic order of the iron sublattice. However, the most 
striking observation is the appearance of two extra CF 
peaks at 9.72 ± 0.05 meV and 11.8 ± 0.1 (labelled B and 
C, respectively) meV in the superconducting compound. 
There are no structural or magnetic phase transitions 
and the amount of PrOF impurity phase is far too small 
to explain them. The temperature dependence of their 
intensities confirm that all the transitions represent tran- 
sitions from the ground state. However, the intensity of 
peak A decreases much faster than the intensities of B 
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FIG. 2: (Color) Inelastic neutron scattering spectra of (a) 
PrFeAsO and (b) PrFeAsOo.87Fo. 13 (circles) measured on IN4 
with an incident energy of 17meV and an average elastic 
wavevector, Q, of 0.85 A -1 . The non-magnetic background 
is approximately given by the spectra of LaFeAsO (trian- 
gles). The solid lines in the panels (a) and (b) are fits to 
Gaussian lineshapes convolved with the instrumental resolu- 
tion, (a) In PrFeAsO, the CF excitations from the ground 
state are centered at 3.58 ± 0.01 meV at 1.5 K. The inset 
shows the temperature evolution of the CF excitation ener- 
gies, (b) In PrFeAsO0.87F0.13, A, B and C label CF excita- 
tions from the ground state, with energies of 2.78 ±0.01 meV, 
9.72±0.05 meV, and 11.8±0.1meV, respectively, measured at 
1.5 K. The inset shows the probabilities as a function of dop- 
ing of the rare earth site having zero (solid line), one (dashed 
line), or two (dotted line) fluorine ions as nearest neighbors on 
the oxygen sublattice assuming a random distribution. The 
solid circle is the ratio of the intensities of the ~ 3meV CF 
excitation in PrFeAsO0.87F0.13 and PrFeAsO. 
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FIG. 3: Temperature dependence of the Van Vleck suscepti- 
bilities derived from the inelastic neutron scattering data for 
the crystal field transitions at ~ 3.5 meV (A) and ~ 9.7 meV 
(B). The temperature dependence of the intensity of peak C 
is similar to peak B. 

and C. As seen in Fig. 2(b), between 1.5 K and 50 K, the 
intensity of peak A falls by a factor of four whereas the 
intensities of peaks B and C remain almost the same. 

In order to quantify this observation, we have fit the 
measured data with a set of three Gaussian peaks and 
converted them to static susceptibilities in absolute units, 
using the Kramers-Kronig relations. The results of the 
temperature dependence of these fits are presented in Fig. 
3. The susceptibilities of both the A and B transitions 
show typical Van Vleck behavior but their different tem- 
perature dependencies unambiguously indicate that these 
transitions belongs to two different rare-earth sites with 
different CF potential or charge environments. 

The wide separation between the peaks at A and B 
show that fluorine doping strongly affects the local elec- 
trostatic potential felt by the rare earth ions. Given that 
there are four oxygen nearest neighbors to the rare earth, 
up to five different crystal field spectra are possible due 
to configurations ranging from zero to four fluorine near- 
est neighbors. Assuming a random distribution of fluo- 
rine atoms at this doping level would result in 57% of 
the rare earth sites with no fluorine nearest neighbors. 
Additionally 34% of the rare earth sites would have one 
fluorine nearest neighbor. The remaining 9% of the rare 
earth sites would have two, three, or four fluorine near- 
est neighbors with decreasing probability. Any tendency 
of the fluorine ions to cluster would alter these ratios. 
In particular, it would increase the fraction of rare earth 
sites with no fluorine ions and redistribute the remain- 
ing ratios. In fact, a comparison of the spectral weight 
of the peak at ~ 3.5 meV between the superconducting 
compound and the parent compound yields a ratio of 
58 ± 4%, which is consistent with a random distribution 
fluorine dopants. Although it should be confirmed by 
measuring this ratio as a function of x, this value is suffi- 
ciently precise to be strong evidence against any substan- 
tial clustering of fluorine ions on the oxygen sublattice. 

The widths of the CF transitions in both the par- 
ent and superconducting compounds are quite small and 



have a Gaussian line shape. The ~ 3meV peak width 
(FWHM) at 1.5 K is 0.61 ± 0.02 meV in PrFeAsO and 
1.54 ±0.03 me V in PrFeAsO0.s7F0.13, after correction for 
the instrumental resolution. This indicates that there 
is some broadening from chemical disorder, probably 
produced by lattice strains due to longer-range fluctu- 
ations in the configuration of fluorine neighbors. The 
fact that we are seeing well defined, sharp CF transi- 
tions in the superconducting compound is evidence that 
the distribution of structural or electronic defects around 
the praseodymium sites is small and that we are dealing 
with two well-defined charge environments. 

Another example where fluorine doping pro- 
duces multiple rare earth environments in a 
flFeAsOx-aFamaterial is for R =Nd. The INS 
data for NdFeAsOo.85Fo.i5 measur ed on the Merlin 
spectrometer at 7 K and 50 K and at initial energies 
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FIG. 4: (Color) (a) Inelastic neutron scattering spectra of 
NdFeAsOo.gsFo.ismeasured on Merlin with incident energies 
of 15meV (circles) and 35meV (diamonds) at an average Q 
of lA _1 and temperatures of 7K (blue symbols) and 50 K 
(red symbols). Open symbols are non-magnetic scattering 
from LaFeAsO0.85F0.15. In the panels (b)-(e) the inelastic 
neutron scattering data measured on the IN6 with an inci- 
dent energy of 3.1 meV and an average elastic Q of 0.65 A -1 
for NdFeAsOo.8 5 Fo.i5(b), (d) and PrFeAsO . 8 rF . 13(c), (e) at 
10 K (b), (c) and 30 K (d), (e). The dotted line in (b), (d) and 
solid line in (c), (e) is the elastic nuclear scattering. In (b) and 
(d), the dashed line is the quasielastic Gaussian line shape fit 
and the solid line is the sum of elastic nuclear, magnetic and 
small linear background contributions. 
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15meV and 30meV are shown in Fig. 4(a). For Nd 3+ 
ions in a point group symmetry lower than cubic, i.e. 
tetragonal, orthorhombic, or monoclinic, the J = 9/2 
ground state multiplct breaks up into five Kramers 
doublets, so there can be at most four CF transitions 
from the ground state. There are clearly five peaks due 
to ground state CF transitions and thus, while it is 
difficult to assign each peak to a particular rare earth 
site, the data are consistent with the picture presented 
above, in which the local electrostatic potential of the 
rare earth ion is modified by fluorine doping. 

These observations bear some similarity to recent NQR 
measurements showing the existence of two charge envi- 
ronments in underdoped i?FeAsOi_ x F x for R =La and 
Sm[5T]. 75 As NQR is sensitive to the the electric field 
gradient, i.e., the second-degree CF potential, acting on 
the arsenic sites, which is also affected by changes in the 
local fluorine ion distribution. However, Lang et al ar- 
gue that the two spectral components they observe in 
the superconducting phase is due to local electronic or- 
der on the iron layer, because of the x-dependence of the 
spectral weights. It is more difficult to model the effect of 
different fluorine ion configurations because of the greater 
distance of the arsenic ions from the oxygen/fluorine sub- 
lattice and the longer-range of the second-degree CF po- 
tential, at least in ionic environments. It is clear that 
their data is more affected by inhomogeneous broaden- 
ing than the neutron data making the spectral weight 
ratios more difficult to determine, so we cannot rule out 
their interpretation, but it seems unlikely that two such 
similar probes should both find evidence of two charge 
environments with completely different origins. 

The final issue we wish to comment on is the effect 
that rare earth substitution has on superconductivity. As 
stated earlier, the fact that the rare earth substituents 
are isovalent to lanthanum mean that the increase in 
T c is cither due to chemical pressure from their smaller 
ionic size or due to the 4/ magnetic moments coupling 
to the iron magnetism and influencing the pairing in- 
teraction. The high resolution neutron scattering data 
shown in Fig 4(b-e) show that the low energy magnetic 
fluctuations are very different in PrFeAsO0.s7F0.13 and 
NdFeAsO0.85F0.15- In the case of neodymium, there 
is strong quasielastic scattering indicating a magnetic 
ground state as expected for a system of Kramer's dou- 
blets. On the other hand, the praseodymium sample ex- 
hibits no magnetic signal in this energy window, indicat- 
ing unambiguously a nonmagnetic singlet ground state. 
Therefore, even though the superconducting transition 
temperature is nearly the same for both materials at op- 
timal doping, the rare earth ground states are very differ- 
ent in both materials. This makes it unlikely that the 4/ 
magnetic moment is involved in the superconductivity in 
these materials. Thus, we conclude that it is the effect 
of chemical pressure due to the lanthanide contraction 
that is responsible for enhancing T c when lanthanum is 



replaced by another rare earth element. 

In conclusion, we have measured the crystal field ex- 
citations in i^FeAsOi^aTa; using inelastic neutron scat- 
tering and established the existence of two charge envi- 
ronments for rare earth sites in nearly optimally-doped 
PrFeAsOi_ x F x and NdFeAsOi-^F^compounds that are 
due to a random distribution of fluorine ions, rather than 
electronic phase separation on the iron layers as pro- 
posed in an earlier NQR study [2"T]. Measurements of 
low-energy magnetic fluctuations reveal that Pr 3+ and 
Nd 3+ 4/ elecrons have nonmagnetic and magnetic CF 
ground states, respectively, from which we infer that the 
4/ magnetic moments are not responsible for the nearly 
identical enhancement of superconductivity in these com- 
pounds compared to LaFeAsOi-^Fa, and conclude that 
a more likely candidate is chemical pressure produced by 
the lanthanide contraction. 
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